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We present a heuristic derivation of the first passage time exponent for the integral of a random 
walk [Y. G. Sinai, Theor. Math. Phys. 90, 219 (1992)]. Building on this derivation, we construct 
an estimation scheme to understand the first passage time exponent for the integral of the integral 
of a random walk, which is numerically observed to be 0.220 ± 0.001. We discuss the implications 
of this estimation scheme for the nth integral of a random walk. For completeness, we also address 
the n — <x case. Finally, we explore an application of these processes to an extended, elastic object 
being pulled through a random potential by a uniform applied force. In so doing, we demonstrate a 
time reparameterization freedom in the Langevin equation that maps nonlinear stochastic processes 
into linear ones. 
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I. INTRODUCTION 



We investigate the general random walk x{t) obeying the equation of motion 

d"-x{t) 



= v{t), (1) 



where rj is white noise with zero mean and unit variance and x{0) = Xo- 

Let us begin with n = 1. The first passage time is the time it takes for the walk to reach zero. When the walk has 
no bias, as above, there is no definite time to expect such an event and the distribution is a power law. To find the 
first passage time distribution p{t)dt, start an ensemble of random walkers at Xo > and at time t — 0. Whenever a 
random walk reaches zero, it is removed from the ensemble. Let Pt{x) be the number density of walks at time t and 
positive X. Pt{x) is a solution of the diffusion equation with absorbing boundary condition, -Pt(O) = 0. More precisely, 

Pt(a;) = — =(e 2t _ g 2* ) 



«2a;„-^(^e-^) (2) 
dx V27rt 

at long times. 

Let g(t) be the integral of this probability distribution over positive values of x. This is the probability that a 
random walk at time t has not crossed zero. The first passage time distribution is then given by 

p{t)dt ^ -^dt = . (3) 

^ ' dt ./2^ t ^ ' 



Although we will not show it, this result is universal for all symmetric walks. Any random walk which is equally likely 
to move forward as backward by a given amount has a first passage time distribution with the same asymptotic form, 
j . This is the Sparre- Anderson theorem [1]. Note that x"^ sets the time scale. 
' The previous discussion is one of the few first passage time problems where an exact solution may be easily found by 
^ . solving a Fokker-Planck equation with absorbing boundaries. Extensions of this method to more complicated random 
■ - - ' walks, such as the second-order random walk described by 



<^x 
'dt 



2 - Vit), (4) 



exist [2]. However, we investigate the first passage time distribution for this walk, and for any n, in a different way. 
At long times, the first passage time distribution for these processes is a power law, 

1 dt 

(5) 
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As computed above, /3i = i. /32 is known to be j [3]. The others arc most likely not exact fractions and are nontrivial 
critical exponents of certain statistical models. Wc first present numerical results for n =^ 2,3,4,5. Then, after 
presenting a heuristic derivation of /32, we make a quantitative estimate for the shift in first passage exponent for 
n > 3. Our analytical results will draw heavily from the theory of Levy flights [4]. We also address separately the 
n = 00 hmit [5]. 

Finally, within zero temperature mean field theory, we demonstrate that the first passage time exponent in the 
n = 2 case is the avalanche size exponent for a dynamically elastic extended object, like a crack front or interface, 
being pulled through a random medium by a uniform applied force at a special point in the parameter space [6]. 
Generically, for these nonequilibrium systems there is a transition from an overall stationary phase to an overall 
moving phase as the applied force is increased toward a critical value. Increasing the applied force from the static 
side triggers local motion of the interface for some finite amount of time as long as the applied force is held fixed from 
the time the toppling starts until the time the toppling stops. The distribution of amount the interface has moved 
during an "avalanche" event, i.e. the avalanche size, gives us information about the continuity (or discontinuity) of 
the depinning transition and, therefore, has been a focus of study over the years [7]. 



II. LEVY FLIGHTS 



A Levy variable Qi is a random variable with a power law distribution for large q, 

P{q)dq oc^^, (6) 

with < (3 < 2. The variance < > is infinite for these distributions and /? is called the Levy exponent of q. 

A Levy flight is a random walk with each step a Levy variable. It is the sum of many independc;ntly drawn Levy 
variables. Because of the infinite variance, a Levy flight is an irregular walk dominated by the few largest jumps [4]. 

Let LffiQ) be the distribution of Q = J2i=i 1i- The Fourier transform Liq{k) is the iVth power of the Fourier 
transform of P{q). If P{k) were twice differentiable at zero, then ijv(fc) ~ (1 - bk"^ + Oik^))^ « 6"''^'= for b real 
and positive. This is the central limit theorem, and it applies when the second moment of a distribution is finite. For 
a Levy variable, the second moment of the distribution is infinite, and P{k) has a cusp at zero for < /3 < 2. The 
form of the cusp may be determined as follows: for small k, the Fourier transform is approximately the integral of 
the distribution over the first wavelength A = or 

When there is a cusp, the limiting distribution is not a Gaussian, but has the following form [4] 

Zjv(fc) =e-*^l*=l'' fc>0 (8) 



By rescaling k, b can be made into a pure phase. For the case where P{q) is symmetric about zero, the Fourier 
transform is real and b — I. It is clear that this distribution has tlic^ correct form near fc = 0, and increasing is 
equivalent to a rescaling of k. The limiting distribution gets wider without changing shape — so it is a fixed point of 
convolution. 

We will need one special nonsymmetric distribution, the distribution of a flight composed of only positive steps. In 
this case, Ln[Q) is zero for Q negative, therefore L]s[{k) must be analytic in the lower half-plane. In other words, 
L{k) for fc < is the analytic continuation of L{k) for A; > 0. As fc is rotated to —k through negative imaginary 
values, k^ acquires a phase e~^'^^. We conclude that 

b* = be-'""^ (9) 



b = e'^f^/^ 



For the special case /? = 5 , b= ^ . 
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And finally, if we were to compute the first passage time distribution for the process, 



!=«*). (10) 

where £,{t) is Levy noise (a Levy variable symmetrically distributed about zero), then we would find a first passage 
time exponent of |. Even though the Levy flight is irregular, the position does not cross zero any faster than it does 
in the nearest-neighbor random walk. This is the Sparre- Anderson theorem once more, which we refer the reader to 
Ref. [1] for the details. 



III. NUMERICAL RESULTS 



To obtain the first passage time exponent /3„ numerically is not as easy as it might appear. Direct numerical 
integration of the equation of motion becomes more cumbersome with increasing n. To efficiently simulate the 
equation, we have calculated the frcc-spacc propagation kernels for n — 2,3, 4,5 (Appendix A). In this context, the 
propagation kernel is the probability distribution of x{to + tg), given initial values x{to), x'{to), ■■■ , a;^"'~^^(io) and 
final values x'{to + ts), ... x^'^~^^to + ts). We first generate the highest-order time derivative, x^"~^^(to + tg), and 
then new values for each lower-order time derivative until x(to + tg) is updated. 

The next time step tg is chosen so that the variance for the next value of x will be smaller. The ratio of the new 
variance to the old is x- X should be small so that the walk does not become negative then positive within one time 
step. The probability of this occurring is exponentially small in the inverse of x- With this algorithm, if x{to) is large, 
the time steps are large as well. The number of time steps required for the simulation only grows logarithmically as 
a function of the first passage time. Near the end of the simulation, when x is small, the time steps become small as 
well, and accuracy is not sacrificed. 

The first passage time distributions for n = 2,3,4,5 are shown in Fig. 1 on a log-log scale. Each plot contains 
approximately 10^ samples in bins of doubling size. Table I contains the linear regression values of the exponents. 
The results are independent of x over the range [0.025, 0.005], indicating that x is small enough. 



IV. THE FIRST PASSAGE TIME EXPONENT FOR n = 2 



For the n = 2 case, x{t) is the integral of a random walk. In other words, the variable that executes the random 
walk is x'{t), the velocity. To exploit this fact, we introduce a new time counter i and divide the time axis into 
intervals Ati between the points where the velocity crosses zero. Then the interval sizes Ati are first passage times 
for an ordinary random walk. 

There is one complication. Referring back to Eq. (3), the time scale until a zero-crossing is the initial value squared. 
So right after a zero- crossing, the initial value is zero. Since the probability distribution is singular at zero, the random 
walk reaches zero again instantly, and then again, infinitely many times. This is a well known property of random 
walks — they jiggle about every value before moving on. Since this affects only the distribution of the smallest times, 
we cut off the distribution of times near zero by imagining the system on a lattice. Now we have a finite and discrete 
set of time intervals Ati with each interval distributed with the Levy exponent /3i = ^. These time intervals are 
independently distributed because the velocity undergoes a simple random walk. 

To compute the first passage time for the position, we must compute the total time t = 5^^]^ Ati until the position 
becomes negative. N is the first zero-crossing of the velocity which happens at a negative value of the position so we 
are actually slightly overestimating t by summing N intervals. We can control for this by summing N — 1 intervals 
instead, which would be a slight underestimate. This estimate of t, be it over or under, should not affect the first 
passage time exponent since we do not expect the last position step (corresponding to the partial time step) to be 
arbitrarily large. 

To continue we need to know two things — 1) the distribution of N, f{N), and 2) given N, the probability 
distribution of t. 

The first problem is actually no problem at all. Since the velocity is just as likely to go up at the start of a 
time interval as it is to go down, the position is as likely to be positive as negative. Therefore, the position is a 
symmetric walk on the zero-crossings (actually, a symmetric Levy flight with independent steps distributed with Levy 
exponent 1/3 using scaling arguments) . By the Sparre- Anderson theorem [1] , the distribution of the first-passage- A'' 
is asymptotically the same as for an unbiased, simple random walk, and 

m) ~ (11) 
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for large N. 

And now for the second problem. The distribution of the sum of N independent At's, Lisf{t)dt, is a Levy distribution. 
The fact that the steps are distributed with Levy exponent /3i = ^ and that all of the At's are positive fixes the 
Fourier transform of the Levy distribution. The inverse Fourier transform can be computed exactly in this case. For 
large t, 

^ y ^ , N dt 

LN{t)dt^ -^e-^j. (12) 

This distribution is small for A''^ > t, and it may be approximated by for < t. Using scaling, analogous 

results may be derived for /3 ^ 5, where a closed form does not exist. 

Now, the time intervals we arc adding are not really independent. This is because we are restricting attention to 
Levy flights that end on the iVth step and this an atypical sample of all Levy flights with N steps. However, the 
properties we actually use from the previous distribution — namely that it is small for N'^ > t and the appropriate 
power law for large t — are the scaling laws for the sum of almost any collection of N Levy variables, where N is large. 
So, although we do not prove it, we assume that these properties hold for the correct LN(t)dt, or the distribution of 
the time elapsed for those Levy flights that end on the Nth step. This is why we refer to our derivation as heuristic. 

All the ingredients are now in place for computing /32- There is a probability f{N) for any value of N, and given A'', 
we know the probability distribution of t. So to find the total distribution of t, we sum up the conditional distributions 
L^it) weighted by the probability of N, or 

00 

P2{t)dt=Y,fWLN{t)dt. (13) 
iV=0 

Approximating the sum with an integral and using the approximate scaling form for Lff{t), 

P2{t)dt^ [ dN f{N)LN{t)dt^ ^^dN. (14) 
Jo Jo N2 t2 

We obtain the following asymptotic form, 

P2{t)dt = ^^, (15) 

so that P2 = J- This result is unchanged when A^ is shifted by one unit, so that the last step is of no consequence for 
the exponent as there are typically many zero- crossings, making our derivation self- consistent. 
Y. G. Sinai has rigorously computed /J2 [3]; we have used some similar methods. 



V. THE FIRST PASSAGE TIME EXPONENT FOR n > 3 



The next random walk we consider is the random surge process, governed by 

d^x da , . 

The acceleration is now an ordinary random walk. After a time t, the acceleration, velocity and position scale as 
o-typ ^ ^'typ ^ t^i and aJtyp ~ t^/"^, respectively. We proceed as we did for the n = 2 case. Once again, we divide 
the time axis by the zero-crossings of the velocity. The time intervals Atj are now distributed with the Levy exponent 
02 = \- The sum of N independent Af's is approximately zero for t < A'^, and for larger t, is approximated by 

T , . , N dt 

LN{t)dt^—-. (17) 
ti t 

Proceeding casually, we may think that, as before, the distribution f{N) has Levy exponent of ^ since that result 
is universal for all symmetric walks. We perform an integral analogous to Eq. (14) and obtain a first passage time 
Levy exponent /Js = |. However, our numerical simulations yield a Levy exponent of 0.220 ± 0.001. The exponent is 
closer to 7 than g. Clearly, there is something wrong. 
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The method fails because there are now correlations between the Axi steps, where Axi is the change in x during 
a Ati step. In the previous n = 2 case, the velocity undergoes a simple random walk and so its sign is random after 
a zero-crossing. In this n = 3 case, the velocity is the integral of a random walk. It is a difTerentiable function and 
so almost certainly changes sign when it hits zero. Consequently, the Axi steps alternate in sign. In addition, as the 
acceleration drifts about, the Ax, steps increase in size. The acceleration at the beginning of each next time interval 
is larger on average than the previous one, making the next Ati and, therefore, the next Axi larger than the last. 

A large number of alternating, increasing steps will reach zero quickly. From the numerical simulation, with each 
oscillation there is a definite probability of reaching zero which is almost constant. In other words, the distribution 
f{N) is not a power law as in the n = 2 case, but an exponential, 

f{N) ~ e-^=*(2iv+i) ^j^j^ ^ Q ^ Q Q2. (18) 

See Table II and the corresponding Fig. 2. The number of velocity zero-crossings must be odd because the position 

can become negative only after the velocity has become negative. With this exponential distribution, the number 
of time intervals one must add up is actually quite small. However, adding a finite number of independently drawn 
Levy variables distributed with Levy exponent j produces a flight with Levy exponent | , so it is quite surprising that 
adding N of them together, where has a finite mean, leads to anything other than an exponent of j. 

And yet, adding a small number of correlated Levy variables does shift the exponent. To quantitatively describe the 
correlation in magnitude, we need to determine how the initial acceleration sets the time scale for each Ati interval. 
We saw from Eq. (3) for the n = 1 random walk that the time scale is the initial value squared. Therefore, the time 
scale between two acceleration zero-crossings is the initial acceleration squared. Because the time between velocity 
zero-crossings is comprised of many acceleration zero-crossings, the time scale between velocity zero-crossings is fixed 
in the same way. 

This argument is somewhat general, so it is nice to verify that it is correct. The coefficient in the formula for 
the first passage time distribution was computed in Rcf. [2] for the case of the n = 2 random walk by solving the 
Fokker-Planck equation with the appropriate boundary conditions. 

P2{t\Vo,Xo^ Q)dtr~. ^— = { — )2— (19) 

for large t. Note that t^yp ~ Since the acceleration in the n = 3 case is exactly the same as the velocity in the 
n = 2 case, we can translate this result into n = 3 language by replacing v with a. Therefore in the n = 3 case, the 
square of the initial acceleration determines the time scale until the next zero-crossing. 

So the square of the initial acceleration at each time step determines the scale of the next time step. Since the 
acceleration is an ordinary random walk, its square has a size proportional to the total elapsed time. Therefore, 
At = T X q where g is a unit Levy variable and T is the total elapsed time. To describe this correlated process, we 
set the units of time so that Ati is a unit size Levy variable qi. The next time step, At2, is no longer imit sized, but 
has a magnitude determined by the square of the acceleration, or, equivalently, the total elapsed time. So to find At2 
we take Ati and multiply it by separate, independent, unit Levy variable q2- To find Ats, we take the total elapsed 
time, Ati + ^t2, which is the expected square acceleration, and multiply by q^. In equations, 

Ati = qi (20) 

At2 = {Ati)q2 

At3 = {Ati + At2)q3 

Ati = {Ati + At2 + Atz)q4. . 
The total time t that has elapsed after N steps is then 

N 

t = qi^{l + qi). (21) 
i=2 

The asymptotic distribution is extracted when qi » 1. This means wc need to determine the distribution of 
products of Levy variables. Taking the natural logarithm, the product becomes a sum. Define Z2 = ln{At2), 
fi = ln{qi), and /2 = ln{q2)- We use the following approximate distribution for both Levy variables. 
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W,= ^^f (22) 
P{q)dq = q <Co. 

with Co = 1. Keep in mind, in order to compute fis, the Levy exponent a of the qi variables equals j. In terms of the 
transformed variables, the distribution of /i and /2 is given by 

P{f) = ae-"f . (23) 

If qi and q2 had different Levy exponents, ai > 0:2, the distribution of the sum would be 

r dhP{h)P{z2 - h) = "1"^ (e-"^ - e-«^ ^^) . (24) 
Jo [0:2 -ai) 

For large Z2, P{z2) ~ Q!2 e~"=^=, i.e. the smaller exponent dominates the tail. However, when ai = a2 = a, we have 

P(z2) = a^^ae-""^ . (25) 
In terms of qi and 52, the probability distribution of gig2 is [8] 

P(AWAi2)=a^^^. (26) 

In general, the probability distribution for i independent Levy variables multiplied together is 

^(^^) = (T^^r'e-"^* . (27) 
The (i — 1)! is required for normalization. Transforming back to the original Ati variables, 

P(At,)^(AtO=(-^^-^^. (28) 

We substitute i = N and a = ^ and this distribution becomes Li^{t) — the probability distribution of the iVth 
passage time of the velocity. Unlike the iVth passage time of the acceleration, L m (t) is not a Levy stable distribution 
because of the correlations. More specifically, LN{t) is the distribution of a product of N Levy variables YiQi- The 
expression for the A'^th passage time we derived earlier is n(9i + 1)) which differs in subleading terms from the previous 
expression. 

Note that L^lt) coincides with the A^tli passage time distribution for the position in the n = 2 case. Both the 
velocity for n = 3 and the position for n = 2 are integrals of an n = 1 random walk. Observe that the distribution 
1^2 (i) is not the same as the distribution of the sum of two variables distributed as Li{t) — there is an additional 
logarithm. We determined L2{t) numerically and it agrees with this form (Fig. 3). For higher N, there are N — 1 
factors of ln{t). 

Before estimating we must make one last argument in the form of an approximation. Up until this point, 

LN{t)dt assumes that each step is free to be larger than the previous one. This is because we have not yet taken into 
account the conditional dependence of the correlation on the position being positive for the walk to persist. In fact, 
the generic situation when the steps are ever-increasing and alternating in sign is for there to be only one step. This 
means that walks that survive more than one step are special. For these walks, the second step is an exceptional pick 
from the distribution. The easiest way for the walk to survive is for the second step to be approximately the same 
size; as the first step. In other words, the size of the first step exits off the distribution of the second step. The result is 
that the second step shares the distribution of the first. For the third step, the position is positive, and now its time 
interval is distributed the same as the second step would have been. Generalizing our argument, any two complete 
time intervals can simply be lumped together as one and we conclude that the distribution of the total after N steps 
is the same as the distribution of the step sizes after ^ steps. In other words, LN{t)dt is not the correct distribution 
for A'' steps, but LN/2{t)dt is. We present evidence for the even-numbered complete steps being typically small in 
Fig. 4. For n = 3, the second step is, on average, about five times smaller than the first time interval. For n = 4, the 
second time interval is typically twice as small as the first. 

So we now have all the ingredients to estimate the Levy exponent P^. We arrive at 
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Ill 

ps{t)dt ~ J2 ^j^{ln{t)re-^^^''dt (29) 

JV=0 

= e"^^^e^'"«dt, (30) 

which gives /Jf^* = ;j — ^-4— • Substituting the numerical value A3 = 0.79, we find the Levy index 0.199. This result 

is a shift in the right direction from our original estimate, but it is still too large. 

Figure 4 suggests that our estimation scheme may work better for n > 3 and so we extend the method to n = 4, 5, .... 
As before, we define the two Levy flights of the time and position intervals between velocity zero-crossings. For n > 3, 
the size correlation factor for the time interval Levy flight turns out to be the same factor of the elapsed time as in the 
n = 3 case with each step retaining only the memory of the highest-order derivative. To be more precise, the initial 
acceleration at the start of each time interval for the n = 3 case is replaced with the highest-order random variable. 
We ignore any other memory effects. Then, for the nth-order process, the step-size distribution for the Atj intervals 
are governed by the passage time exponent for the (n — l)th-order process. Referring back to Eq. (28), i = N and 
(X = Pn-i- In addition, the correlation in the sign of the Ax^ still persists for the nth-order random walk so the 
distribution f{N) remains exponential with a decay constant that we determine numerically. Therefore, the nth Levy 
exponent will be a small perturbation about the (n — l)th first passage time exponent with a shift of — /?„_ie~^^". 
For n = 4 and n = 5, we observe that our estimation scheme yields more accurate results. See Table IL 

As in the n = 2 case, for n > 3 we consider the last partial time step to be a full one. Given that L]^/2it)dt is the 
more appropriate distribution, the bulk of the first passage time for a given is taken up during the last two steps 
(up to subleading corrections). So the these last partial time steps should be distributed just as the first passage time. 
We have numerically verified this to be the case. 

Majumdar et. al. [5] have constructed a completely different approximation scheme to compute the first passage 
time exponent for these higher-order random walks. While we refer the reader to Ref. [5] for the details, their scheme 
is rather accurate for (^2 and /Js. 



VI. THE n = 00 LIMIT 



As n increases, we also expect the decay constant A„ to grow. In fact, in the limit that n becomes large, the 
first passage time exponent reaches a limiting value. We consider the solution to Eq. (1) with initial conditions 
a;(0) = x'(0) = ... = a;"(0) = 0. It is 

X{t) = [ dtn [ dtn-1... I dtiTjih), (31) 



which is a linear function of r] and therefore is the inner product of r] with a kernel. In this case, 

;-t')^"" 



c{t) = [ 
Jo 



dt'^^TI^^Vin (32) 



It is easy to see that this expression satisfies Eq. (1). 

When solving for the time when x{t) first becomes zero, it is permitted to rescale x{t) by any finite number, even if 
that number is a function of time. This does not change the set of points where x{t) is zero. Rescaling by so that 
y{t) = and dropping constant factors, 

y{t) = jl^dt'il-jr-'v{t'). (33) 

In the large n limit, (1 — y)" is indistinguishable from an exponential in the region where it has the most weight 
{t' < nt). The result is that 

y{t) = - f dt'e^riit') . (34) 

Majumdar it et al. [5] have noted that the n = 00 random walk has the same first passage time exponent as a very 
different system — the solution of the two-dimensional diffusion equation with a random initial condition. 
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(35) 



p{x,y,0) = 'r]{x,y) . 
In polar coordinates, the value of p(0, t) is given by 

P{0, ^) = ^ / ^'^e- - -^vir, 0) , (36) 

where we have used S(x)5{y) = r^^6{r)d{9). The angular integration leaves the integral in terms of a new random 
variable, 77' (r) = J dOri{r, 9). We now have 

p{Q,t) = j dr^e-'^^riir) (37) 

Transforming variables to zt = , we arrive at 

1 1 



^ y due-^r]{u), (38) 



where we used the fact that d{r) = S{u)2^/u so that rj' [r) = ri{u)\^ui . 

Comparing this kernel with the kernel for the n = 00 random walk, wo sec that they arc indeed the same. As a 
result, the distribution of times for the diffusive field to change sign at one point is the same as the distribution of 
the first passage time for the n = 00 random walk for long times, where the initial conditions are irrelevant. 



VII. PHYSICAL SYSTEMS 



Consider a planar crack front moving through rough, brittle material. In the quasistatic regime, the planar crack 
front is a line with long-range static elasticity. There is a uniform stress applied to the material, which drives the 
crack front forward. But there are also local, random variations in the fracture toughness which the crack front needs 
to overcome to move forward. These randomly varying forces arc pinning forces. The competition between the global 
pulling and local pinning forces as mediated by the elasticity of the crack front determines its dynamics. For small 
values of the applied force, the pinning forces eventually dominate and the crack front remains stuck. For some finite 
value of the applied force, the pinning forces are not strong enough to keep the crack front in place and it moves 
forward. The boundary between these two phases is the depinning transition. 

Approaching the depinning transition from the static side, a small increase in the applied force causes a tiny 
portion of the crack front to move forward. This portion causes a few more regions to move forward as well, only 
to be eventually stopped by more strongly pinned regions. The only role of the applied force here is to induce the 
initial motion of the tiny portion. It remains fixed until the local motion of the crack front stops. Near the depinning 
transition, there exists a sequence of these discrete, localized avalanche events with distribution of sizes r. The 
avalanche size is defined as the total amount of the crack front that has moved forward during these events. A power 
law distribution of avalanche sizes as the force approaches the critical force indicates a continuous, second-order-like 
depinning transition. 

Near the depinning transition, the motion of the crack front is jerky. We therefore model its dynamics in terms of 
discrete space and time. Within an infinite-range model, where each crack front segment is coupled equally to every 
other, the spatial degrees of freedom along the crack front arc averaged out. We define yt to be the total number 
of segments of the crack front that hop forward at discrete time t. Given that we are only concerned with the few 
segments that are on the verge of hopping forward, yt is Poissonian distributed. In this infinite-range limit, and in 
the dissipative regime, the mean of yi is determined by the number of segments that have hopped at the previous 
time step only. In addition, the mean of yt fluctuates statistically until it reaches zero, whereupon the motion stops. 
In the limit of large avalanches where things are changing slowly with time, the equation of motion for the avalanches 
at the depinning transition is 

% = y-^n{t) (39) 

in the Ito representation (the time derivatives are forward differences). When y(t) reaches zero, the avalanche is over. 
We refer the reader to Ref. [6] for a more detailed derivation. 
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So far, wc have neglected any effects hke elastic waves that are indeed present along the crack front. If a piece 
of the crack front moves forward, that motion creates an extra transient force on its neighbors as the elastic wave 
propagates along the crack front. The extra transient stress is called a stress overshoot. This effect is different than 
pure inertia, which would be an extra transient stress on the segment itself after it jumps. However, we demonstrate 
that the two effects are similar in this infinite-range model for a particular value of the stress overshoot [6]. If we 
were to take into account elastic waves along the crack front, the mean of yt depends on several previous time steps. 
In other words, there are higher-order corrections to Eq. (39) that are usually irrelevant in the long time limit. The 
most significant terms are the first and second time derivatives, and the first time derivative is the only relevant term. 
If we fine-tune the coefficient of ^ to zero, wc arrive at a tricritical point. Here, the second-order derivative is the 
most relevant term and we obtain the equation of motion 

^ = yiv{t). (40) 

The variable we are interested in is 

rtf 



Jo 



dty{t), (41) 



the avalanche size. We emphasize that tf is determined by the boundary condition y{tf) = 0. 

In order to find the distribution of t, we use the size of the avalanche as a new path-dependent time coordinate. We 
transform time on each particular path in the path- integral differently, in a way that depends on the history of y{t). 
This transformation alters the density of ?/(t)'s at any given time, mixing together 2/(t)'s from different times so that 
they appear simultaneous. Because of this property, it is difficult to perform the transformation in the Fokker-Planck 
equation; but the transformation is easy in the path integral. 

A reparameterization of time, even a path dependent one, does not affect the answers to questions that do not 
involve the time explicitly. The probability for a stochastic walker to arrive at point B from point A is independent of 
the global time, but the probability for a stochastic walker to be at point B at time t is not. Similarly, the avalanche 
size of a given walk does not involve the time lapse of the walk as the upper cutoff of the integral is determined by 
the boundary condition on y{t) and not the time explicitly. 

Choosing the new time coordinate r to tick at the new rate. 



and using the property of Gaussian noise, ?7(t(<)) = the (Ito) equation pair 



^ = y~'v{t) (43) 



dy 

— = V 
dt 

become 

^ = *) (44) 

dy 

2 

Replacing y with x ~ ^^ completes the transformation. The zeros of y are also the zeros of x so the value of t when 
X is zero is then the avalanche size. We therefore find the distribution of avalanche sizes to be the distribution of first 
passage times t for the n = 2 random walk. It has a Levy exponent j. 

We may perform the same transformation on the Fokker-Planck equation directly, although the motivation then 
becomes obscure. Consider the following probability-conserving equation, 

dP dP 1 d^P dP 
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for a new quantity Pt{x,v,T). If we integrate Pt{x,v,T) over all r, wc recover the original Fokker-Planck equation. 
On the other hand, if integrate Pt{x,v,T) with respect to t from zero to infinity, we obtain 

dK _ Id^ vdK_ 
dr 2 dv'^ y dy ' 

where K{v,y,T) = dt Pt{v,y,T). Defining before, gives 

dK _ 15^ _ dl£ 
dr 2 dv"^ dx ' 



(47) 



which is the Fokker-Planck equation for the n = 2 random walk. This demonstrates the equivalence of the two 
problems. 

It should now be clear that the avalanche size exponent for Eq. (39) is ^ since the same; time reparameterization 
may be performed in that case also. This relates the avalanche size exponent to the first passage time of the n = 1 
random walk. 



VIII. CONCLUSIONS 



We have presented a relation between the shift in the first passage time exponent and the decay rate of the 
probability of N velocity zero-crossings for the nth random walk. The method can be exact for n = 2, but not for 

n > 3 because of correlations. More work is needed to estimate the decay constant of f{N) since this would determine 
the convergence rate of the exponents to the n = oo value. In addition, we do not have any bounds on the error of the 
first passage time exponent as the approximation is uncontrolled. We also do not know which types of correlations 
for Levy flights lead to shifts in exponents and which do not. A classification of correlations in terms of these two 
categories might be useful. 

We have opted to slice up time in terms of velocity zero-crossings, instead of using the global time in the Fokker- 
Planck equation. This approach allows us to analyze the higher-order random walks in terms of a one-dimensional 
Levy process since the phase variables are the position and the time only. The remaining variables are subsumed in 
the [n — l)th-ordcr random walk that the velocity undergoes between the zero-crossings. 

With Eq. (42) we have also reparameterized the global time, but in a different way. We use a path-dependent 
time transformation to analyze the avalanche statistics of a nonlinear second-order random walk in terms of the first 
passage time exponent of the linear second-order random walk. This allows us to give a novel physical interpretation 
of the first passage time exponent for the n = 2 case. 

Since the avalanche size does not depend on the global time, it has been known to us that one can find the stationary 
solution of the Fokker-Planck equation corresponding to the nonlinear random walk in Eq. (40) and then arrive at 
the avalanche size exponent after invoking a simple scaling argument [9] . Since this method gives a new derivation of 
the avalanche size exponent, it indirectly gives a new derivation of the first passage time exponent for the n = 2 case. 
This insight may be useful to construct a derivation of /3„ for n > 3. 

Finally, these higher-order random walks could be relevant for more general forms of dynamic stress transfer along 
a crack front, or other extended clastic objects. For instance, one can approach a fine-tuned critical point where the 
coefficients of the first and second time derivative terms are zero. Then the third-order time derivative is the most 
relevant term, and we arrive at a higher-order critical point. While the time reparameterization given by Eq. (42) 
docs not allow us to equate the avalanche size exponent with the first passage time exponent in this case, we remain 
optimistic that either exponent may be relevant for various physical systems. 
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APPENDIX: THE KERNEL OF HIGHER-ORDER RANDOM WALKS 



We begin with the partition function for all higher-order random walks, 

Z = J D[ri{t)]e~^-!o^"^'^'^\ (Al) 

We change variables from r](t) to x{t), where = r]{t), and use forward difference time derivatives. We change 
variables with the insertion D[x{t)]5{^^-^^ — r){t)) to arrive at 

/ D[x{t)]e~^io . (A2) 

Jx{0)=Xo 

We then minimize the action to find the classical equation of motion with vanishing variations on the boundary for 
all derivatives up to the (n — l)th- derivative. In the n = 2 case we find the equation of motion, 

We then impose the following constraints on the x{t) trajectory, a;(0) = Xq, x{T) = Xf, x'{0) = Vg, and x'{T) = Vf, 
yielding 

(3.x /■ - 3.i:„ - r,T- 2r„r) . 

X{t) =Xo+ Vot + ^ ^j;^ 

^ {vf + Vq - §.{Xf - Xq)) ^ ^^4) 

The path integral is quadratic, so the partition function is proportional to the exponential of the classical action. 
Substituting the classical solution into the action and integrating gives the propagation kernel, 

Pt{Vf,Xf\Vo,Xo) = —p^e 2* 

TT t 

e ^ . (A5) 

Wc normalize the propagation kernel to have integral unity. Notice that the kernel is the product of a Gaussian in v 
and a gaussian in x. Given that x{t) is the time integral of v{t) and that the sum of Gaussian variables is Gaussian, if 
the velocity is Gaussian, so is the position. Not only does the position spread faster than the velocity, but it also drifts 
with time with a coefficient that depends on the velocity. This method works for higher n giving us the following 
propagation kernels for n = 3, 4, 5: 

Pt{af,Vf,Xf; ao, Vo, Xo) - le-^(«^-«°)%-^r«''^-""')-^(«^+«"')*)' 

t2 



-^i{x,-Xo)-i{vo+Vf)t+^{af-ao)t^r 



for n = 3, 



1 (Sf-s^)2 

Pt{sf,af,Vf,Xf;So,ao,Vo,Xo) ~ -g-e 



e 



for n = 4 with = s{t), and 



dt 

Pt{uf,Sf,af,Vf,Xf;Uo,So,ao,Vo,Xo) ~ -^e" e-^(^*-f-'''^-i^"^+"°'*^' 

t 2 

.S^((af-ao)-i{sf+So)t+^iuf-Uo)t^)\-^{{vf-Vo)-i{af+ao)t+^{sf-So)t^-Ti^(uf+Uo)t^f 
p-^^^T^{{=Of-Xo)-k{v,+Vo)t+^{af-ao)t^-M'f+So)t^+Tiso(^f-»o)t^f 



(A8) 
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for n = 5 with = u{t). 
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TABLE 1. Tlio first i)a.s.sagt' tiiiu' t'xptnioiits oljtaiui'd from tits to tiic data sliowu in Eig. 1. 


n 




2 
3 
4 
5 


0.005 0.250±0.001 
0.005 0.220±0.001 
0.005 0.210±0.001 
0.005 0.204±0.001 


A. 


TABLE II. The first passage time exponent estimates for the n > 3, which uses the numerical exponential decay constant 
I and the numerical results for Pn-i- This estimate is to be compared with the last column in Table I. 


n 


X A„ 


3 
4 
5 


0.005 0.79±0.02 0.199 
0.005 1.55±0.02 0.210 
0.005 1.92±0.02 0.205 
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FIG. 1. Log-log plot of the probability of a first passage time t (in arbitrary units) occurring within the interval [t/y/2, V^t), 
where the probability has been multiplied by a different constant for each curve so that they do not overlap. The open circles 

denote numerical data for n — 2; the open squares for n = 3; the open diamonds for n = 4; the open triangles for n = 5. The 
solid lines represent the results of the linear regression. The size of the symbols is larger than the error bars. As we only used 
double precision in our simulations, there is an upper cutoff in the first passage time of approximately 10^®. 
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FIG. 2. Log-linear plot of the distribution of number of velocity zero-crossings f{N), where f(N) has been multiplied by a 
different constant for each curve so that they shift and do not overlap. The symbols here are the same as in Fig. 1. For the 
n = 3 curve, we ignore the first data point in the linear regression as there is some memory of the initial conditions for this 
point. 
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FIG. 3. Log-log plot of Li (t) and {t) for the n = 2 case, i.e. the probability of a first (or second) position passage time t 
(in arbitrary units) occurring within the interval [t/\/2, \/2t). The open circles denote numerical data for Li{t); the inverted 
triangles for L2{t). While the solid line represents the results of the linear regression; the dot-dashed line represents Eq. (28) 
up to a proportionality constant, with i = 2 and Ati = t. As we must keep track of even larger absolute values of the position 

for N = 2, the double precision constraint cuts off the tail of the distribution even more so for N = 2 than for N = 1. Once 
again, the symbols are larger than the error bars and there are approximately 10*^ samples. 
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FIG. 4. Plot of the probability of the ratio of the second complete time interval to the first complete time interval, denoted 
as s, occurring within the interval [s, s + 0.01) for the n — 3 (the open squares) and rt = 4 (the open diamonds) random 
walk. Both probabilities eventually fall off exponentially as opposed to being distributed as L2{At2)d{At2), indicating that the 
second step is indeed an exceptional pick from the distribution for the first time interval. Each curve contains approximately 
10® samples. 
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